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Abstract 



At critical coupling, the interactions of Ginzburg-Landau vortices are determined by the 
metric on the moduli space of static solutions. Here, a formula for the asymptotic metric 
for two well separated vortices is obtained, which depends on a modified Bessel function. 
A straightforward extension gives the metric for N vortices. The asymptotic metric is 
also shown to follow from a physical model, where each vortex is treated as a point-like 
particle carrying a scalar charge and a magnetic dipole moment of the same magnitude. 
The geodesic motion of two well separated vortices is investigated, and the asymptotic 
dependence of the scattering angle on the impact parameter is determined. Formulae 
for the asymptotic Ricci and scalar curvatures of the iV-vortex moduli space are also 
obtained. 
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1 Introduction 



Critically coupled Ginzburg-Landau vortices and BPS monopoles are two of the most studied 
examples of topological solitons in field theory [Hj. Vortices are particle-like solutions of the 
abelian Higgs theory in two dimensions and BPS monopoles are solutions of a Yang-Mills-Higgs 
theory in three dimensions. 

At critical coupling (separating the Type I and Type II regimes of superconductivity), 
vortices exert no static forces on each other, and there are static multi-vortex solutions. These 
satisfy the planar Bogomolny equations jH] 

D l( j) + iD 2 (j) = (1.1) 
B + i(H 2 -l) = 0, (1.2) 

and the boundary condition \<f>\ — > 1 as |x| — > oo. Here, Di<p = di<p + iAi<p is the covariant 
derivative of the complex scalar field 0, and B = diA 2 — d 2 Ai is the magnetic field in the plane. 
Taubes showed that an N- vortex solution is uniquely determined by specifying N points where 
is zero OH]. The moduli space of N- vortex solutions, A4n, is therefore the configuration 
space of N unordered points in the plane, which is a smooth 2iV-dimensional manifold. Using 
the first Bogomolny equation, the gauge potential can be eliminated, and the second equation 
can then be written in terms of the gauge invariant field h = log |0| 2 as 

N 

V 2 h-e h + l = 47r^5(x-y r ). (1.3) 

r=l 

The points {y r : 1 < r < N} are the locations of the vortices, where vanishes and h has a 
logarithmic singularity. The boundary condition is h — * as |x| — > oo. Taubes showed that h 
approaches exponentially fast [Hj. 

Static BPS monopoles are solutions in M 3 of the Bogomolny equations 

Bi = D&, (1.4) 

where D$ is the covariant derivative of the adjoint Higgs field and Bi is the Yang-Mills 
magnetic field. For fields of finite energy there is a well-defined monopole number N, and the 
moduli space of iV-monopole solutions for gauge group SU(2) is a AN- dimensional smooth 
manifold. It is not so simple as in the vortex case to say precisely what the moduli signify 
without introducing some additional structures (e.g. Donaldson's rational map), but for well 
separated monopoles there are four moduli associated with each of them. Three specify the 
location of the monopole, and the fourth is an internal phase angle. 

We are interested not just in static solutions but also in time-dependent ones. We sup- 
pose the complete Lagrangian, both for vortices and for monopoles, is the Lorentz invariant 
extension of the static energy function, with a kinetic term quadratic in the time derivatives 
of the fields. In the vortex case, the Lagrangian density is that of the abelian Higgs theory at 
critical coupling 

£ = -\f»»F^ + \d^>D^ - ^(|0| 2 - l) 2 , (1.5) 
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where F^ v = d^Ay — d^A^ and D^ifi = <9 M + iA^ip. The vortices can move at arbitrary- 
speeds less than the speed of light. For monopoles the situation is similar, but the theory is 
non-abelian. High speed collisions of either vortices or monopoles are complicated, involving 
substantial energy transfer to radiation modes of the fields, and amenable only to numerical 
simulation. However, collisions at slow speeds can be treated adiabatically, using the geodesic 
approximation The idea here is that the moduli space of static solutions acquires a natural 
metric by restricting the kinetic terms of the field theory Lagrangian to motion tangent to the 
moduli space, and it can be shown that the geodesic trajectories on moduli space accurately 
model the field theory dynamics of the solitons. An argument for this was given in It 
was justified rigorously for vortex motion and for monopole motion [18J by Stuart. It is 
fairly clear now that the geodesic approximation is the formal non-relativistic limit of the field 
theory dynamics of the solitons, where radiation is neglected. 

Having recognized the importance of the metric on moduli space, it becomes desirable to 
calculate it. This is not so easy. For monopoles there is an explicit understanding of the 
metric only for one monopole, where the metric is flat, and for two monopoles, where it was 
calculated by Atiyah and Hitchin pQ. For vortices there is a general formula due to Samols 
[T3] . which we shall use below. Again, for one vortex the metric is flat, but even for two 
vortices the metric is not known explicitly. 

However, it is possible to calculate the explicit asymptotic metric on the moduli space 
for N well separated monopoles. There are now two approaches to this calculation. The 
first is physically motivated, and not quite rigorous [101 d- The monopoles are treated as 
point-like objects, carrying a magnetic charge and a scalar charge of equal magnitude. These 
charges are regarded as sources for auxiliary linear fields through which the monopoles interact. 
(Monopoles, in reality, are smooth but nonlinear, with a core radius of order 1, but provided 
their separation is much greater than 1, a linearization of the fields appears justified.) For 
monopoles at rest, the magnetic forces exactly cancel the scalar forces, so there is no net 
force. For monopoles in relative motion, the magnetic and scalar forces are not identical, 
because of their different Lorentz transformation properties, and there are net forces which 
cause the monopoles to scatter non-trivially. The effects can be encapsulated in an iV-particle 
Lagrangian with a purely kinetic term, quadratic in velocities. The coefficient matrix of this 
quadratic form defines the metric on the moduli space for well separated monopoles. An 
alternative approach is due to Bielawski, who calculated rigorously the asymptotic form of the 
Nahm data associated with well separated monopoles, and from this calculated the asymptotic 
metric on moduli space [2]. These approaches give the same result, and they are consistent 
with the Atiyah-Hitchin metric for two monopoles, whose asymptotic form was derived in 
The asymptotic iV-monopole metric, like the true metric on the iV-monopole moduli space, 
is hyperkahler, but unlike the true metric it has singularities when the monopoles come close 
together. 

These results on monopoles motivated the present work. Here we give an explicit expression 
for the metric on the iV-vortex moduli space, for N well separated vortices. Furthermore we 
calculate it in two ways. Our first approach is to take Samols' general formula and evaluate 
the quantities occurring there by a method of matched asymptotic expansions. Essentially, 
we solve eq. ()1.3|) for two well separated vortices, calculating the effect of one vortex on the 
other at linear order, and from the solution determine the asymptotic 2-vortex metric. It is 
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straightforward to generalize the 2-vortex metric to the iV-vortex metric. We need to assume 
that for well separated vortices the field h far from the vortex cores obeys the linearization of 
(jl.3j) . namely the Helmholtz equation 

\7 2 h-h = 0, (1.6) 

and that the relevant solution is a linear superposition of the solutions due to the N vor- 
tices separately. Corrections due to the nonlinear terms neglected in eq. 1)1.6)1 are of higher 
exponential order in the separations. However, a careful treatment of this point is lacking, 
and would require a considerable refinement of Taubes' estimate of the exponential decay of 
solutions. 

Our second approach is the more physical. It is a variant of the calculation involving 
point-like monopoles, and regards well separated vortices as point-like sources interacting via 
auxiliary linear fields. A study of the static forces between vortices that are close to critical 
coupling shows that well separated vortices can be regarded as particles each carrying a scalar 
charge and a magnetic dipole moment (thought of as perpendicular to the plane which the 
vortex inhabits) jTHj. For critically coupled vortices, the magnitudes of the scalar charge and 
the dipole are the same, and the static forces due to them cancel. For vortices in motion, 
the scalar and magnetic forces do not exactly cancel, but result in velocity-dependent forces. 
Again there is an effective Lagrangian for two well separated vortices which is purely kinetic, 
and from this the metric can be read off. The extension to iV vortices is as before. 

The asymptotic metric for vortices, which involves the Bessel function K , has some simi- 
larities to the true metric. It has the same isometries, and like the true metric, it is Kahler. As 
for monopoles, the asymptotic metric becomes singular as the vortices approach one another 
closely, since it is not positive definite if the minimum vortex separation is below a certain 
critical value (2.21 to two decimal places in the case N = 2). Of course, the asymptotic metric 
is not valid in this region. 

It remains open to rigorously prove that our formula gives the asymptotic metric on the 
iV-vortex moduli space, but the known results for monopoles make this conjecture plausible. 
Using our formula we can calculate the scattering of two vortices that do not approach close 
to each other. The leading exponentially small expression for the scattering angle can be 
obtained exactly. 

This paper is organized as follows. In Section 2 we obtain the asymptotic 2-vortex metric, 
and its generalization to iV vortices, along with the asymptotic Ricci and scalar curvatures. In 
Section 3 we rederive the metric using the model of vortices as point-like particles. In Section 
4 we discuss the scattering of two vortices using the asymptotic metric. 

2 Well Separated Vortices — Field and Metric 

The key to the metric on A4]y, the N- vortex moduli space, is the equation (J1.3)) . whose 
solutions determine the static iV-vortex fields. It is sometimes convenient to use a complex 
coordinate z = x + iy for a general point in the plane, and to denote the vortex locations 
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correspondingly by {Z r : 1 < r < N}. Eq. f)1.3j) becomes 



A? 



V 2 h-e h + l =4irJ2 5 ( z - Z r) (2-1) 

r=l 

where V 2 = 4~|^=. Around Z r , the function /i(z, z) has the local expansion 

h = log | ^ — Z r \ 2 + a r + -b r (z — Z r ) + -6 r (5 — Z r ) 

+ c r (z-Z r ) 2 -^(z-Z r )(z-Z r )+c r (z-Z r ) 2 + ... , (2.2) 

where a r is real, and 5 r , c r complex. Taubes proved that this series, with the logarithmic term 
removed, is a convergent Taylor expansion. The logarithmic term and the coefficient \ are 
determined by the equation locally, but the remaining coefficients are not. They depend on 
the positions of the other vortices, but not in an explicitly known way. Most important for us 
is the coefficient b r . 

Samols' formula for the metric on M.^ is 

N 

I ill) \ 

(2.3) 



g = 7T ^ ( S rs + 2 ^-j dZ r dZ s 

r,s=l ^ ' ' 



The functions b r obey the symmetry relation 

db , db r 



dZ r dZ v 



(2.4) 



and from this it follows that the metric is not only real, but also Kahler. Invariance of 
the metric under a translation of all the vortices implies that ^2 b r = ^3], arid rotational 
invariance implies that Yl Z r b r is real |13j . 

For well separated vortices, we assume that h is exponentially small except in a core region 
with radius of order 1 around each vortex, and there h has an approximate, local circular 
symmetry. It follows that if the minimum separation of any pair is L 3> 1, then the S rs term 
dominates the metric, and the correction is of order e~ L . The metric is therefore approximately 
flat. 

Let us now concentrate on two vortices, and denote their positions by 

Z x = Z + ae ie , Z 2 = Z - ae ie . (2.5) 

It follows from the symmetry of the 2- vortex field around the centre of mass Z, or from the 
properties of the functions b r mentioned above, that in this case b\ = b(a)e %e and 6 2 = —b\, 
where b(a) is a real function. Samols' formula implies that the moduli space metric is 

g = 2ndZdZ + 7](cr)(d(T 2 + a 2 d9 2 ) (2.6) 

where 

"< ff > = a ( 1 + ^K>))- (2 ' 7) 
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The relative motion of two vortices takes place on the reduced moduli space where Z is fixed. 
This is a surface of revolution. The range of the coordinates is < o < oo and — | < 9 < |, 
with 9 = — | and 9 — | identified. The range of 9 is it and not 2it because the vortices are 
identical. Therefore, the surface is asymptotically conical, rather than planar. 

So far, our exposition has been a summary of known results, but now we show how to 
calculate the leading asymptotic correction to the conical metric. We return to the equation 
(jl.3|) . and consider first the circularly symmetric solution ho for a single vortex at the origin. 
In terms of polar coordinates (p, <p), the equation satisfied by h , for p > 0, is 

d 2 hn 1 dhn . 

-TT- + -- 7^- e + 1 = 0- 2 - 8 
dp z p dp 

The boundary conditions are ho ~ 21ogp for small p, and h — > as p — > oo. The Taylor 

expansion of h — 2 log p about p = involves only even powers of p. For large p, eq. (|2.8|) 
has the linearized form 

d 2 h 1 dh , , , 

the modified Bessel equation of zeroth order, so 

ho(j>) ~ -K (p) , (2.10) 

7T 

where g is a constant. The corrections to this asymptotic approximation are expected to be 
suppressed by order e~ p . By numerical integration of the nonlinear equation ()2.8|) . it has been 
determined that q = —10.6 [TSJ. Recently, Tong has given an argument involving dualities 
in string theory which strongly suggests that q = — 27r84, in agreement with the numerical 
result [201 ■ There is as yet no direct proof of this using (|2.8|) . 

Next, let us consider the perturbation of the solution h due to other, distant vortices, still 
assuming that one vortex, which we label as vortex 1, is precisely at the origin. Let us write 
h = h Q + hi, where hi is small in the neighbourhood of vortex 1. The linearization of eq. (jl.3|) 
implies that 

(V 2 -e ho ) h = 0. (2.11) 

The operator acting on hi has no singularity at the origin, so hi is smooth there, and the 
logarithmic singularity of h is carried entirely by ho. Since ho is circularly symmetric, we can 
separate variables and write 

h(p, <p) = ho(j>) + h(p,ip) 

j oo 

■=ho{p) + -fo(p) +^2[f n (p) cos n<p + g n (p) sin n<p) , (2.12) 

" n=l 

where /„ obeys the equation 

d 2 fn + l_dfn_( eka + <S) fn = 0i (2M) 



dp 2 p dp \ p 2 

and g n obeys the same equation. f n is nonsingular at p = and has a series expansion 
fn = o. n p n + . . .. Similarly, g n = (3 n p n + . . .. The expansion (|2.12|) is consistent with the general 
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expansion of h around vortex 1, that is, (|2.2|) with r = 1 and Z\ = 0. By identifying the terms 
linear in p, we find that b\, the coefficient we are interested in, is given by b\ = ax + i(3\. 
From now on, therefore, we just consider equation (|2.13J) for fx, that is, 



d fi 1 dfx ( h 1 \ 

+ -#- e feo + - /i = 0. (2.14) 



<ip 2 p dp \ p 2 



For large p, this simplifies to 



^4+ if -(i+ii/^o, 

ap z p dp \ p^ 



which is valid for 1 C p <C L, where L is the distance from the origin to the next-nearest 
vortex. Note that the difference between the coefficients in eqs. 1)2.14)1 and 1)2.15)1 is e h ° — 1, 
which is smooth, finite, and exponentially localized. We therefore suppose that the asymptotic 
form of the solutions of ()2.14j) are exact solutions of (j2.15j) . This is supported by the results 
used in various examples of scattering theory, and which follow from Levinson's theorem jl]; 
however a result of the precise type we require, involving perturbations of Bessel's equation, 
appears to us not to have been established. Eq. (J2.15)) is the modified Bessel equation of first 
order, whose general solution is a linear combination of the functions Kx(p) and ii(p). 

Let us now assume that there is just one other vortex, vortex 2, whose location (in Cartesian 
coordinates) is (— 2a, 0), with a 1. In this case, h is reflection symmetric under <p t— > —ip, 
so in the Fourier series for hx all the functions g n vanish. In particular, bx — ax, and is real. 

In the region far from both vortex centres, the equation (Jl.H)) linearizes to 

V 2 /i-/i = (2.16) 
and is solved by the linear superposition of the fields due to each vortex separately 

h(p, <p) = ^Ko(p) + -K f v / 4^ + 4apcos<p + p 2 ) . (2.17) 

The argument of the second Kq function is the distance to vortex 2 from the point with polar 
coordinates (p, ip). By separation of variables, the general solution of the Helmholtz equation 
(12. 16)1 . regular at p = and with the reflection symmetry <p i— > — <p, is a linear combination 

of the functions I n (p) cosmp. The function K 4cr 2 + 4crp cos <p + p 2 j is such a solution 

(whereas K (p) is not, being singular at p = 0), so 

oo 

K (V 4 ^ 2 +4apcos<p + p 2 ) = k I (p) + 2 ^ fc rJn(p) cos(rep) (2.18) 

n=l 

for some real constants k n . 

Note that an important special solution of ()2.1fij) is e x = e pcosip , and its expansion 

oo 

= J (p) + 2^J n (p)cos(n<p) (2.19) 

n=l 
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defines the functions I n {p)- Combining the series for the exponential function with trignometric 
identities, one can compute the leading terms in the series expansions for I n . These are also 
given in standard references, e.g. [Hj. It is sufficient for us to record that 

Jo(p) = l + ..., h{ P ) = \p + ... . (2.20) 

We can now return to (J2.18)) and determine k\, the coefficient we need. The Taylor expan- 
sion (in p) of the two sides gives 

K Q {2a) — Ki(2cr) p cos ip + . . . = k + kipcosp> + . . . , (2.21) 

where we have used the identity K\ = —Kq, and the results above for J and I\. So k = K {2a) 
and 

k 1 = -K l {2a). (2.22) 

With this result, we can now match the Fourier expansion of (j2.17|) . the linearized field h 
due to the two vortices, valid outside their cores, with the Fourier expansion of h = ho + hi 
near vortex 1. In the range 1 <C p <C 2a, we find 

?-(Ko(p)+K (2a)I (p)) - ^K 1 (2<t)I 1 (j>)cxxi<p+... 

7T V / 7T 



±K (p) + i/ (p) + f\(p) cos<p + . . . (2.23) 

7T L 



Therefore, fi(p) has the asymptotic form 



h(p) = -^ Kl (2a)h(p) (2.24) 

7T 

and there is no K-y(p) piece. The further terms on both sides of (|2.23J) involve cosrop with 
n > 1, and could be determined from higher order terms in the Taylor expansion (|2.21|) . but 
we do not need these. 

The last step is to extrapolate the function f\ into the core region of vortex 1. It is rather 
remarkable that this can be done, because the equation satisfied by fi, namely (j2.14j) . is not 
a standard equation, and the coefficient e h ° is not known explicitly. However, one solution of 
(j2HH is known. It is 

This can be verified by differentiating (|2.8|) . the nonlinear equation for ho. The interpretation 
of this solution is that it corresponds to the translational zero mode of vortex 1. If the centre 
of that vortex is infinitesimally translated by e in the x-direction, then the field of vortex 1 
becomes h = ho(p — ecostp) to first order in e. So h = ho + h\ where hi = —efi cos<p and 
fi = Since h = ^K (p) asymptotically, it follows that /i = —&Ki(p) asymptotically. 

Similarly, since h ~ 21ogp for small p, it follows that f\ ~ 2/p for small p. By contrast, the 
solution fi of (j2.14j) that really interests us has the asymptotic behaviour — — Ki(2a)Ii(p) for 
large p, and the finite linear behaviour b\p for small p, where b\ is to be found. 
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Now, equation ()2.14j) has a Wronskian identity 

dfi d fi 

P ^ dp ~~ ^ dp ' = constant ' ^ 2 - 26 ^ 

relating the two solutions /i and f\. Using the asymptotic forms of f\ and fi, and the 
Wronskian identity for the modified Bessel functions 

we deduce that the constant in (j2.26|) is ^-i^i(2cr). Evaluating ()2.26|) near p = we deduce, 
finally, that 

&! = ^K x {2a) . (2.28) 

We can now use this result to calculate the 2-vortex metric. In the above calculation, 
vortex 1 was at the origin and vortex 2 at (— 2a, 0). From ()2.5|) we see that Z = —o and 
9 = 0, so 

b(a) = J^i(2<r) . (2.29) 
Therefore the prefactor 77 in the metric (|2.fi|) is 

77(a) = 2tt fl - ^#o(2<7)^ , (2.30) 

where we have used ()2.7|) and the identity K[(s) + K\(s)/s = —Kq(s). The complete asymp- 
totic 2-vortex metric is 

g = 2vr dZdZ + 2tt (\ - ^K (2a)^ (da 2 + a 2 d# 2 ) . (2.31) 

We shall investigate the geodesies of this metric in section 4, and hence determine how vortices 
scatter. 

To extend ()2.31|) to the asymptotic iV-vortex metric is not hard. Let us use the complex 
coordinates of the vortices Z r and introduce the notation Z rs := Z r — Z s . The flat part of the 
metric ()2.3|) can be reexpressed as 

N 

71 dZ r dZ r = NndZdZ + — ^ dZ rs dZ rs , (2.32) 

r=l r^s 

where Z is the centre of mass coordinate 

Z = ^(Z 1 + Z 2 + ... + Z N ). (2.33) 
Note that the differentials dZ rs are not all linearly independent. 
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To find the remaining part of the metric, we need to find b s and its derivatives. The 
solution of the Helmholtz equation ()2.16|) becomes a linear superposition of the fields due to 
the N vortices. The asymptotic matching of h in the neighbourhood of the s-th vortex can be 
carried out as before. This leads to the following expression for b s that is a linear superposition 
of the effects of the other N — 1 vortices, 

<f V- TS nrr ,N Z s 



6 ' = ^£*i(i z »i)TfV ( 2 - 34 ) 



J sr 



Each term is the obvious generalization of ([2.28)1 . combined with the orientational phase factor 
Z sr /\Z sr \ which reduces to e l6 for two vortices. 
Because of translational invariance, 

£lr°. < 2 - 35 > 

r=l 

SO 

db s ^ db s 

dz s = ~^dz: r (2 - 36) 

(no summation over s). For r ^ s, we find, differentiating (|2.34|) with respect to Z r and 
keeping Z r fixed, that 

H = £W-D- P.37) 
Eq. (|2.37|) combined with (|2.36|) gives 

|^ dZ r dZ s = £^J2 K o(\ Z sr\) (dZ r - dZ s ) dZ s . (2.38) 

r,s=l 7 r^ts 

Since K (\Z sr \) = K (\Z rs \), we symmetrize over the contributions of these two terms, obtain- 
ing 

q£ dZ r dZ s = K o(\Z rs \) dZ rs dZ rs . (2.39) 

r,s=l r r^s 

Putting these ingredients together, we obtain our final expression for the asymptotic N- 
vortex metric 

g = N7idZdZ + 7iJ2(^-^M\Zrs\)) dZ rs dZ rs . (2.40) 

For two vortices, located at the points (|2.5)L this reduces to (j2.31|) . Since the coefficients in 
the asymptotic metric depend only on the magnitudes of the vortex separations, it is clear 
that it is translationally and rotationally symmetric. The structure of the metric as a small 
perturbation of the flat Euclidean metric becomes clearer if we eliminate the centre of mass 
coordinate Z using (|2.33|) : 

q 2 

g = n^2 dZ r dZ r - — ^ K (\Z rs \) dZ rs dZ rs . (2.41) 

r r^s 
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One way to see that this metric is Kahler is to note that eq. (|2.4|) is satisfied, since ()2.37j) 
and ()2.36|) imply that ^jh is real and symmetric. More explicitly, the asymptotic Kahler form 
is 

u = % J^L d Z AdZ+^J2(^- £l K °(\ Z rs\)) dZ rs A dZ rs . (2.42) 

r^s 

Since the 1-forms dZ, dZ rs are closed, one finds that 
ia 2 K'(\Z \) 

duo = — — > ° VI r . sW (Z rs dZ rs A dZ rs A dZ rs + Z rs dZ rs A dZ rs A dZ rs ) = , (2.43) 
167T \Z rs \ 

so uj is closed. The Kahler potential is 

N 2 

vr V Z T Z T - q - V K Q ( I Z rs I ) . (2 .44) 

z — <• n £ < 

r=l r^s 

The Kahler form is of direct interest in certain non-relativistic models of vortex dynamics 
[T3*] . Such models have first order dynamics in time, and it is conjectured that slow vortex 
dynamics is well approximated by a Hamiltonian flow on the iV-vortex moduli space, where 
the symplectic structure is precisely this Kahler form. Clearly the closure of uj is crucial for 
this to make sense. 

The curvature properties of soliton moduli spaces are of some interest. For example, the 
scalar curvature of M.n is relevant to quantum iV-soliton dynamics ^2], while in the case 
of monopoles, Ricci flatness of M.^ was the key property exploited in Atiyah and Hitchin's 
construction of the metric for N = 2. In order to compute the asymptotic Ricci tensor for the 
N- vortex metric (j2.41j) . it is convenient to write g as 

g = ^ g rs dZ r dZ s = ^ n(8 rs + h rs )dZ r dZ s , (2.45) 

r,s r,s 

and work up to linear order in the perturbation h. It is a standard result in Kahler geometry 
[2T] that the Ricci tensor associated with g is 

r,s ' b 

where G is the determinant of the hermitian coefficient matrix g rs . In this case, 

G = det7r(I + /i) =n N (l + tih + • ■ •) 
log G = N log 7r + h rr H 

r 

= iVlog7r-^VA' (|Z r -Z s |) (2.47) 
from (|2.4ip . Equations ()2.46|) and (|2.47|) together with Bessel's equation imply that 

R=£lJ2 K ^ - Z s\)dZ rs dZ r8 . (2.48) 

r^s 
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One sees that the iV-vortex Ricci tensor is asymptotically positive semi-definite, its two- 
dimensional null space being tangent to the translation orbits in M.^ (that is, a vector is null 
if and only if it generates a rigid translation of the N- vortex system). Tracing R, one obtains 
the scalar curvature, 

Seal = Y^^Rrs = £sJ2 K ^ Z r ~ Z ^ + ' ' ' ' ( 2 - 49 ) 

r,s r^s 

whence one sees that M.n is asymptotically scalar positive. It is an interesting open question 
whether the true metric on A^jy has similar curvature positivity properties. The numerical 
results of Samols for N = 2 suggest that it may |14j . 

3 The Point Source Formalism 

In this section we rederive the asymptotic 2- vortex metric from a more physical viewpoint. 
The idea is that, viewed from afar, a static vortex looks like a solution of a linear field 
theory with a point source at the vortex centre. We will see that the appropriate point 
source is a composite scalar monopole and magnetic dipole in a Klein-Gordon/Proca theory. 
If physics is to be model independent, the forces between vortices should approach those 
between the corresponding point particles in the linear theory as their separation grows. This 
idea, which originated in the context of monopole dynamics ^U], has already been successfully 
used to obtain an asymptotic formula for static intervortex forces away from critical coupling 
[To] . The present application is somewhat more subtle since we are required to analyze the 
interaction between point sources moving along arbitrary trajectories. We handle the problem 
perturbatively: using a mixture of Lorentz invariance and conservation properties we obtain 
expressions for a moving point source and the Klein-Gordon/Proca field it induces correct up 
to acceleration terms. From these we construct the interaction Lagrangian for one moving 
point source interacting with the field induced by another. This Lagrangian is purely kinetic, 
i.e. quadratic in velocities, and hence may naturally be reinterpreted as the energy associated 
with geodesic flow on the asymptotic 2-vortex moduli space. The extension to iV-vortex 
dynamics is entirely trivial. 

In this section 2 ) denotes a space-time point. x° = t is the time and 

x = (a; 1 ,^ 2 ) denotes a spatial point. To linearize the abelian Higgs theory (jl.5j) . we choose 
the gauge so that the scalar field is real. Since vortices have nontrivial winding at infinity, 
this requires a gauge transformation which is singular at the vortex centre. This need not 
concern us since we seek only to replicate the local, far field behaviour of the vortex in the 
linear theory. In this gauge, the vacuum is = 1, so we define <fi = 1 + ip and linearize in ip. 
The resulting Lagrangian density is 

C = l -d^d^ - - \f, u F^ + l -A,A» + m/> - fA, (3.1) 

where k is the scalar charge density and j the electromagnetic current density. These will be 
chosen to replicate the vortex asymptotics. The corresponding field equations are 

+ = k (3.2) 
+ = f (3.3) 
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where □ = d^d^ = d\ — V 2 , and we assume that j is a conserved current, d^j^ = 0. 

In this gauge, the scalar field of a single vortex located at the origin is <fi = exp(~/i ) where 
ho satisfies (|2.8j) . so 

</> = l + \h + ...~l + ±K {\x\) (3.4) 

Z Z7T 

for large |x| by eq. (j2.10|) . Hence we seek a point source k so that the solution of ()3.2j) is ip = 
^K (\yi\). Since the static Klein-Gordon equation (Helmholtz equation) in two dimensions 
has Green's function K , 

(_V 2 + l)K (|x|) = 27r5(x), (3.5) 

one sees that 

K = q5(x) . (3.6) 

For a static vortex, the time component of the gauge potential A$ vanishes. The asymptotic 
behaviour of its spatial components Ai is determined by the first Bogomolny equation (jl.lj) 
which, on linearization, implies 

d 1 i/;-A 2 + i(d 2 ip + A 1 ) = 0. (3.7) 

Hence 

A = (A 1 , A 2 ) = (d 2 , -di)1> = -fkx Vif (|x|) (3.8) 

Z7T 

where we have introduced k, the unit vector in a fictitious x 3 -direction orthogonal to the 
physical plane. It follows that the point source which reproduces the asymptotic vortex gauge 
field in eq. (13. 3J) is 

(j°,j) = (0,-gkx V5(x)). (3.9) 

The physical interpretation of ()3.6|) and ()3.9|) is that the point particle corresponding to a 
single vortex at rest is a composite consisting of a scalar monopole of charge q and a magnetic 
dipole of moment q. We shall refer to this composite as a (static) point vortex. 

The interaction between two arbitrary (possibly time-dependent) composite sources (k(i), j(i)) 
and (k(2), j(2)) in this linear theory is described by the Lagrangian 

Lint = J d 2 x (k ( i)V>(2) - 3(1) A (2)») ( 3 - 10 ) 

where {^u\,Au\) are the fields induced by {Ku^ju)) according to the wave equations (|3.2|) . 
(I3.3|l . This is obtained by extracting the cross terms in J <i 2 x£ where (k, j) = (ft(i), j(i)) + 
( K (2), 3(2)) and (ip, A) = (V'(i), ^4(i)) + (V ; (2), ^.(2)) by linearity. Although (j3.10j) looks asymmetric 
under interchange of sources, it is not, as may be shown using fJ3.2J) . ()3.3|) and integration by 
parts. If the sources are chosen to be static point vortices, that is, translated versions of (|3.fij) 
and ()3.9j) . one finds that L int = 0, so static point vortices exert no net force on one another at 
critical coupling, in agreement with the nonlinear theory. 

We seek to compute L int in the case where the two sources represent point vortices moving 
along arbitrary trajectories in M 2 . To do so, we must construct a time- dependent point source 
representing a vortex moving along some curve y(t), say. The construction is guided by two 
principles: first, in the case of motion at constant velocity, the source should reduce to (|3.6|) . 
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(|3.9j) in the vortex's rest frame; second, for any trajectory the vector source j, which represents 
the vortex's electromagnetic current density, must remain a conserved current. The result will 
be correct up to quadratic order in velocity and linear order in acceleration. 

It is straightforward to calculate Lorentz boosted versions of the sources (jH.fijl . ()3.9j) . Let 
£ M denote the rest frame coordinates and assume that at time t = 0, the point vortex lies at 
the origin, x = 0, and is moving with velocity u. Then, by decomposing x and £ into their 
components parallel and perpendicular to u, one finds that at this time, = y(u)x" (where 



7(u) 



u 



is the usual contraction factor) while £ — x . Hence 



x-u 

y(u) u 



X • u 
X u 

u 



(x • u)u + 



(3.H) 



where the ellipsis denotes discarded terms of order u A or greater. We shall not persist in so 
denoting these terms. Rather, we shall include an ellipsis only where further neglible terms 
have been dropped to obtain the given expression. Since k is a Lorentz scalar, the boosted 
scalar monopole is re(x) = ft s tatic(0 — T° interpret this delta function as a distribution 

on the x-plane, note that for any test function /, 



efx/(x)5(£) = / d 2 £ 



dx 



f(x(Z))6(Z) = (l-~u 2 )f(0) + .... ( 3 - 12 ) 



Therefore 5(£) = (1 — ^u 2 )5(x) and so, at t — 



«(x) = q(l - -n 2 )5(x) . 



(3.13) 



The boosted dipole is more subtle, since j itself transforms as a Lorentz vector. The rest 
frame source is 

Jstatic = (j(°o), j( 0) ) = (0, -qk x V^(^)) . (3.14) 

To obtain the laboratory frame source, one must perform a Lorentz boost on this with velocity 
— u. Explicitly, 



J°(x) 



J"(x) 
J X (x) 



«t(«) i(o)(£( x )) = u - J(o)(£( x )) + 

1 

2 

^)(€(x)). 



7(^ ) (^(x)) = (l + - M 2 )jJI ) (^x)) + 



(3.15) 

(3.16) 
(3.17) 



We may combine ()3.1(j)l and ()3.17j) into a single equation for j by using the same polarization 
trick as in (|3.11j) . yielding 



j(x) = j(o)(^(x)) + -(j(o)(£(x)) • u) u 



(3.18) 



Now 



= (0) V(l - ^n 2 )5(x) 



in 2 )V-iu(u-V) 



5(x) + 



(3.19) 
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Hence 



j°(x) = g(k x u) • V5(x) + . . . (3.20) 
j(x) = -q(l - ^u 2 )k x V5(x) + -g[u(k x u) ■ V+ (k x u) u • V]5(x) 

= -qk x V5(x) + g(k x u) u ■ V5(x) . (3.21) 

By replacing x by x — y(t) and u by y(t) in ()3.13|) . ()3.20|) and (|3.21j) we obtain expressions 
for the instantaneously Lorentz boosted point vortex travelling along an arbitrary trajectory. 
In particular, 

Jboost = g((k x y) • V, -k x V + (k x y) y • V)5(x - y) . (3.22) 

Note that 

dXost = g(k x y) ■ V5(x - y) + . . . ^ (3.23) 

so jboost is not the moving point source we seek. Rather, we must add a correction j acc to jboost 
of order y to enforce current conservation. Such a term will vanish in the case of motion at 
constant velocity and hence does not conflict with the required Lorentz properties of j. We 
make the simplest choice, namely 

Jacc = (0,-gkx y 5(x-y)). (3.24) 

Though we have chosen j° cc = 0, any function of order |y| would do since doj® cc is automatically 
of negligible order. It turns out that this ambiguity has no bearing on our calculation because 
j^. c makes no contribution to L int at order (velocity) 2 . 

To summarize, the point vortex moving along a trajectory y(t) is represented by a com- 
posite point source 

K(t,x) = g(l-i|y| 2 )(5(x- y ) (3.25) 
j(t,x) = ?((kxy)- V,-kx V+(kxy)y- V-kxy)5(x-y). (3.26) 

The second task in the calculation of L int is to construct the fields ("0(2) > ^.(2)) induced by 
the second moving vortex. Were the field equations ()3.2j) and (j3.3J) massless, we could simply 
use retarded potentials, making a suitable expansion in time derivatives. Since they are 
not massless, we need a substitute for this procedure. We handle the problem by introducing 
formal temporal Fourier transforms of the fields and sources, as follows. Let tp(t, x) be the field 
induced by the time varying source n(t,x) according to (|3.2|h and define Fourier transforms 
ijj and k with variable ui dual to t, 



V>(t,x):=/ due iujt ^{u,t) , k(*,x) := / du e iujt lZ{uj,t) . (3.27) 

J~ oo J— oo 

Then ()3.2|) implies 

[-V 2 + (1 -uj 2 )\4) = k, (3.28) 
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so for each value of u, ip{(jJ, •) satisfies the static inhomogeneous Klein-Gordon equation with 
mass yl — uj 2 and source k(uj, •). Comparing with ()3.5|) one sees that (j3.28|) is solved, at least 
formally, by convolution of H(uj, •) with the (suitably scaled) Green's function Kq, 



1 

2^ 



J d 2 x'K (y/l -w 2 |x - x'l) k(u, x' 



(3.29) 



Expanding the Green's function in uj and truncating at order uj 2 yields 



1 

27 



K (|x-x'i)--|x-x'|^(|x-x'|) 



K W, X 



(3.30) 



From this we obtain ■?/> by (|3.27|) . 

lW,x) = i | dV [K (|x-x'|)^,x') + T(|x-x'|)^,x')] 



where we have defined 



T(s) :=-sK' (s). 



(3.31) 



(3.32) 



Note that truncating the expansion in u is, in effect, the same as neglecting higher time 
derivatives of k, and hence eventually of y(t) in our application. No claim of rigour is attached 
to the above Fourier transform manoeuvre. One should regard it as a convenient algebraic 
shorthand for generating a perturbative solution of (J3.2)) . Direct substitution of (J3.31|) into 
(|3.2|) confirms that this really is a solution up to higher time derivative terms (d^K etc.). Since 
equation ()3.3|) is formally identical to ()3.2|) . the vector field induced by a time dependent 
source j is 

A"(t,x) = ±J rfV [K (|x-x'|)^(t,x') + T(|x-x1)^(t,x')] + ... . (3.33) 

Having obtained ("0(2)) A 2 )) induced by the time varying source (ftp), j<2)) we may compute 
the Lagrangian governing its interaction with another source (K,h),jm) by substitution of 
dHISH) and (J3~33j) into (pri"0j) . The result is 



Lint = 7T~ I ^xdVji^odx-x'l) K (1) (t,x)K( 2 )(t,x') - (t, x) J (2 ) M (t, x') 



-Tflx-x'l 



^K(i)(t,x)^K( 2 )(t,x') - d t tf 1) (t,n)d t j( 2 ) t ,(t,x') }(3.34) 



where a total time derivative has been discarded. It remains to substitute the point vortex 
sources (I3.25|) . (J3.26j) for vortices moving along trajectories y(t) and z(t) into (J3.34j) and 
evaluate the integrals. Explicitly, 

Lint = J d 2 xd 2 x[K (1) K i2 ) +j(l) • j(2) - J(°1)J(°2)]^0(|X-X / |) 

J d 2 Kd 2 x![d t K {l) d t K {2) + dtS(i) ■ <9ij (2) - ^'(i)^( 2 )] T (l x - x 1) ( 3 - 35 ) 
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where 



K(i)K( 2 ) 


— « 2 M flvl 

- <? [J- - 2 uy| 


-+- |z| jjo^x — yjo(x — zj 


j(l) • j(2) 


= g 2 [V^-V 2 - 


-(z- V,)(z-V 2 )-(y-V,)(y- V,) 






-z • V y - y • V 2 ]5(x - y)5(x - z 


J(1)J(2) 


= g 2 ((kxy). 


V,)((kxz)-V z )5(x-y)5(x'-z) 


d t K(x)d t K(2) 


= g 2 (y-v y )(z 


■ V,)5(x-y)5(x , -z) 


3tj(l) ' <%j(2) 


= g 2 (y-v,)(z 


• V z )V y - V^(x-y)5(x'-z) 


%'(°1)%'(°2) 


= 0. 





(3.36) 

We have discarded terms of order |y| 3 , |y||z| etc., and systematically used 

V<5(x - y) = -V y S(x - y), V'5(x' - y) = -V^(x' - z) . (3.37) 
All the integrals are thus rendered trivial, so 

Unt = ^{[l-\{\y\ 2 + \i\ 2 ) + V y -V z -{i-V y ){i-V z )-{y-Vy){y-V z ) 
_z • V tf - y • V, - ((k x y) • V„)((k x z) • V z )]K (\y - z|) 
-(y ■ V y )(z ■ V 2 )[l + V, ■ VjT(|y - z|)} . (3.38) 

We now use the following identities, which are easily checked using Bessel's equation (|2.9|) : 

(l + VyV z )K (\y-z\) = (3.39) 
(l + V,-V 2 )T(|y-z|) = -K (\y-z\) (3.40) 

( a .V y )(f3-V z )K (\y-z\) = - a -^l y - z l) - a »^A(|y-z|). (3.41) 

iy z i 

In f!3.41|) . ct, (3 are any vectors with f3 independent of y, and A denotes the function 

AO) := K (s) - 2^^ . (3.42) 
Also, we have used the decomposition of a, /3 relative to the orthonormal frame 

n ii = y ~ z n ± = k x n ii _ ( 3A3 \ 
\y - z l 

Equations (|3.39|) and (J3.4Uj) give useful cancellations in lines 1 and 3 of (J3.38)) . All but two 
of the remaining differential operators (namely, the acceleration terms) are of the form (J3.41|) 
for suitable choice of a and f3. Hence, we may repeatedly apply ()3.4H) . yielding 



L - q2 
Lmt ~ 27 



-i(|y| 2 + |z| 2 ) + (y!l) 2 + (ill) 2 + ^-yllill 



A(|y-z|) + L acc (3.44) 
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where L acc represents the remaining acceleration terms. In fact 

2 2 

Lace = -^(z-V y +y-V z )K (\y-z\) = ^(y-Z).(y-z)K' (\y-z\) 

Z7T ZlX 



2tt 



-,_, !2 M^) + te ,,_ il)2A(|y 



|y - z i 

Substituting (jTOT^i and |H21) into ipCTjl finally yields 



total time , , 
derivative. 



L int = -f\y-z\ 2 K (\y-z\). (3.46) 

Our point particle model of 2-vortex dynamics is completed by adding to L int the usual 
nonrelativistic free Lagrangian for two particles of mass 7r (the vortex rest energy), so 

L = \ (|y| 2 + |z| 2 ) - J^|y - z| 2 K (|y - z|) . (3.47) 

We may define centre of mass and relative coordinates 

R=^(y + z), r=i(y-z) (3.48) 

so that 2 

L = vr|R| 2 + tt f 1 - ^ (2|r|)^ |r| 2 . (3.49) 

Motion under this Lagrangian coincides with geodesic flow on the 2-vortex moduli space with 
respect to the asymptotic metric (j2.31|) of section [21 we simply identify Z = R 1 + iR 2 and 

ae i0 = r l _|_ ^2 

Extension of this treatment to the case of N well separated vortices is trivial due to 
the underlying linearity of our point-particle model. To the standard Lagrangian for N free 
particles of mass ir, moving along trajectories y r (t), r = 1, 2, . . . , N say, one simply adds one 
copy of Li nt , as i n eo ±- Q3-46J) . for each unordered pair of distinct particles, or equivalently, one 
copy of \L mt for each ordered distinct pair. The result is 



\ E i^i 2 - h S K ^ yr ~ yMr ~^ 2 - (3 - 50) 

r r^s 



The asymptotic formula for the iV-vortex metric readily follows. Defining complex coordinates 
Z r = y\ + iy 2 , their differences Z rs = Z r — Z s , and holomorphic 1-forms dZ rs = dZ r — dZ s , one 
sees that 

q 2 

g = n^2 dZ r dZ r - — ^2 K (\Z rs \) dZ rs dZ rs , (3.51) 
which coincides with (12.41)1 . 
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4 Two- Vortex Scattering 



The relative motion of two vortices, in the geodesic approximation, is determined by the purely 
kinetic Lagrangian 

L = ± V (a)(a 2 + a 2 8 2 ) (4.1) 

where the ranges of a and 8 are as in section 2. Samols calculated the function b(a) and hence 
r){o) numerically, and using this found the geodesic motion of two vortices [Hj. The geodesic 
motion has two constants of integration: the energy E, which is L itself, and the angular 
momentum £, which equals r](a)a 2 8. Using these, one can find d6/da, and the scattering 
angle can be determined by integration. It depends on the impact parameter a, given by the 
motion of the vortices as they approach from infinity. There the energy is rev 2 and the angular 
momentum is 2'nav. 

The result for the scattering angle as a function of impact parameter agrees well with 
numerical simulations of 2- vortex scattering using the complete field equations [TT) I15j. The 
motion is repulsive, and the scattering angle increases monotonically as the impact parameter 
decreases, from zero when the impact parameter is infinite, up to | in a head-on collision. 
The field dynamics begins to significantly differ from the geodesic motion only if the vortex 
speeds exceed about half the speed of light. 

Now that we have obtained the asymptotic form of the 2-vortex metric, with r](a) given 
by eq. (|2.30jh it is possible to estimate the asymptotic scattering. Unfortunately, the exact 
relation between scattering angle and impact parameter for this metric is given by a rather 
intractable integral. Rather than evaluate this numerically, we adopt a different strategy. 

Since our asymptotic metric is only valid at large vortex separations, it can only be used 
with confidence to evaluate the scattering angle at large impact parameter. Here the vortex 
trajectories are almost straight, and the scattering angle small. We can find the approximate 
scattering angle by a perturbative calculation. It is convenient to use Cartesian coordinates 
x = a cos 8 and y = <7sin#. The Lagrangian becomes 

L= l -r l (a)(x 2 + y 2 ) (4.2) 

where a = \/ x 2 + y 2 . The equations of motion are 

n'(a) 

r](a)x + -^—^(xx 2 + 2yyx - xy 2 ) = (4.3) 
2a 

r]{a)y+^-^{yy 2 + 2xxy-yx 2 ) = 0. (4.4) 
la 

We may assume that the motion is approximately along the line x — a, with y increasing 
from — oo to oo at approximately constant speed v. (This is in fact the trajectory of one of 
the vortices; the other moves in the opposite direction along the line x = —a.) The initial 
value of x is taken to be zero, and by calculating the small change in x we can calculate the 
scattering angle. We work to leading order in the small quantity exp(— 2a). rj' is of first order 
in exp(— 2a) along the trajectory, so at this order x, is given by the term proportional to y 2 
in eq. (|4.3|) . and the coefficient rj multiplying 'x can be approximated by 2n. rfx is negligible. 
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y is of order exp(— 2a) too, but the consequent change of speed in the ^/-direction along the 
trajectory can be neglected. Thus, it is a sufficiently good approximation to take the solution 
of eq. (|4.4J) to be y = vt, and to simplify eq. (|4.3J) to 



The total change in x is therefore 



47T 

and the scattering angle, assuming it is small, is 



x = °^^. (4.5) 
47r a 



av 2 f°° ri'(a) , , . 

Ax = / dt , (4.6) 



Q = -Ax. (4.7) 

v 



Expressing Ax as an integral over y, using y = vt and a = \J a? + y 2 , and also rf{o) ■ 
^K x (2a), we find 

B = tlf «>ME±n dy (4, 
'-oo ^a 2 + y 2 
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,2 



q 2 d 



2n 2 da 

Remarkably, this gives the simple result 



/oo 
K (2y/a* + y>)dy. (4.9) 
-oo 



a 2 

= ^_exp(-2a). (4.10) 
2tt 

The integral can be understood as follows. The planar Helmholtz equation with source at 
the origin 

(-V 2 + 4)x = 2tt(5(x) (4.11) 



has the exponentially decaying solution x — -^o(2\/a^~+~^). Integrating the equation with 
respect to y we find that \ = K (2^/x 2 ~-\-y 2 ) dy satisfies 

dx 2 



+ 4 )x = 2ir5(x) (4.12) 



and hence K (2^/ ' x 2 + y 2 ) dy = | exp(— 2\x\). 

If we substitute the numerical value of q, we can compare the dependence of scattering 
angle on impact parameter with Samols' result. The agreement is good for a > 2. This is 
shown in Fig. ^ 

We expect corrections to this calculation. These are partly due to the neglected terms in 
(I4.3J1 and (|4.4|) . and partly due to corrections to our asymptotic metric. 
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Figure 1: Scattering angle O against impact parameter a for 2- vortex scattering in the geodesic 
approximation. Dashed curve: Samols' numerical implementation; solid curve: our perturba- 
tive approximation. 
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